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ABSTRACT 



Aims. The basic mechanism responsible for the widespread condensation of MgS in the outflows from carbon rich stars on the tip of 
the AGB is discussed with the aim of developing a condensation model that can be applied in model calculations of dust formation in 
stellar winds. 

Methods. The different possibilities how MgS may be formed in the chemical environment of outflows from carbon stars are explored 
by some thermochemical calculations and by a detailed analysis of the growth kinetics of grains in stellar winds. The optical properties 
of core-mantle grains with a MgS mantle are calculated to demonstrate that such grains reproduce the structure of the observed 30 
yum feature. These considerations are complemented by model calculations of circumstellar dust shells around carbon stars. 
Results. It is argued that MgS is formed via precipitation on silicon carbide grains. This formation mechanism explains some of the 
basic observed features of MgS condensation in dust shells around carbon stars. A weak secondary peak at about 33 . . . 36 yum is 
shown to exist in certain cases if MgS forms a coating on SiC. This new feature seems to have occasionally been observed. 

Key words, circumstellar matter - stars: mass-loss - stars: winds, outflows - stars: AGB and post-AGB 



1. Introduction 

Magnesium sulphide is a wide-spread dust component in dust 
forming stars with carbon rich element mixture. Since Mg and 
S are abundant elements, this dust species forms a major com- 
ponent in the dust mixture injected by dying AGB stars to the 
interstellar medium. Any attempt to model the dust production 
by AGB stars and their contribution to the interstellar dust pop- 
ulation requires a knowledge of the production mechanisms at 
least of the abundant dust species. Since for MgS a condensation 
model for this species is presently lacking, we try to develop in 
this paper such a model. We will explain, why we believe that 
MgS condenses in stellar outflows as a mantle on SiC grains. 

Magnesium sulphide was detected as dust component in cir- 
cumstellar dust shells enshrouding some carbon stars and in the 
environment of two planetary nebulae by Forrest et al. (1981) 
by observing a broad and prominent emission band in the far- 
infrared spectrum which apparently is centred around 30 fim. At 
the instant of its detection the nature of the carrier of this emis- 
sion band was not immediately clear, except that the feature can 
only be a solid state absorption band. Goebel & Moseley ( 1985 ) 
proposed MgS as carrier of this feature. MgS has an infrared 
absorption feature just in the wavelength region of the newly de- 
tected feature and chemical equilibrium condensation calcula- 
tions of Lattimer et al. (1978) for carbon rich element mixtures, 
originally performed for identifying possible supernova conden- 
sates, had shown that MgS should be an abundant solid phase 
at temperatures of the order of 600 K or lower. A companion 
paper to the paper of Goebel & Moseley ( 1985 ) by Nuth et al. 
( 119851 ) showed newly determined optical properties of a number 
of solid phases, MgS, CaS, FeS, SiS 2 , FeS 2 , and Fe 3 C, that could 
be responsible for the 30 fim band and which were compared 
to observations in the paper of Goebel & Moseley. Though the 
match between the observed feature and the absorption band of 
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MgS was not perfect, the identification of MgS as carrier of the 
30 fim band has been widely accepted since that time, though 
other interpretations are occasionally discussed (e.g. Papoular 
120001 Grishko et al. 120011 Volk et al. 120021) . 

In order to allow to include a calculation of the MgS emis- 
sion feature in model calculations of radiative transfer in circum- 
stellar dust shells, and thus to enable a quantitative comparison 
between observations and theory, Begemann et al. (1994) deter- 
mined optical constants of Mg. t Fei_ r S (0.9 > x > 0). They ob- 
tained good agreement between the observed emission band pro- 
file and the calculated spectrum in the 30 /mi region for a radia- 
tive transfer model of the circumstellar dust shell of IRC+10216, 
which left little doubt with respect to the interpretation of the 
30jum emission band as being due to MgS. 

A number of observational and theoretical studies of the 
30 fim emission from dust enshrouded AGB and post-AGB stars, 
and from the environment of planetary nebulae appeared since 
that time (e.g. Omont [T993l Omont et al. |1995t Yamamura et al. 
[19981 Jiang etal. [19991 Szcze rbaeta l. [1999l Hrivnak et al. r2000l 
Hony et al. [20021 V olk et al. 120021 Ho ny & Bouwmann 2004; 
Lagadek et al. 120061 Zijlstra et al. 120061 Leisenring et al. 120081 . 
Good examples for the emission band can be found in Hony et 
al. ( 120021 ), Volk et al. J20021 ) for galactic objects and Lagadek 
et al. (120061 ) and Zijlstra et al. ( 120061 ) for the Magellanic clouds. 
They all show unequivocally that MgS formation is a common 
phenomenon for carbon stars on the tip of the AGB. 

The formation mechanism of this MgS is not clear. It is al- 
ready speculated in the first papers by Nuth et al. (1985) and 
Goebel et al. (1985) that MgS forms via a surface reaction on 
carbon grains and this scenario was favoured over MgS con- 
densation as a separate dust species. A labory investigation on 
MgS condensation was conducted by Kimura et al. (2005 ), but 
this gives no direct insight in the formation process. No further 
discussion of this problem seems to have appeared in the astro- 
physical literature so far. In this paper we will show that MgS can 
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Fig. 1. Stability limits of solid phases formed from the most 
abundant elements in chemical equilibrium for a carbon rich el- 
ement mixture (ec = 1-2 eo). The dashed line shows the P-T- 
trajectory of a stationary dust driven stellar wind with a mass- 
loss rate of M = 10" 5 M Q yr -1 
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Fig. 2. Concentrations in the gas phase of S, Si, and Mg atoms 
and their abundant molecular compounds with H, C, O, and S 
for a cooling sequence with P = 10~ 1() bar if formation of solid 
SiC in chemical equilibrium is allowed for (full lines) and if the 
formation of SiC is suppressed (dashed lines) 



only be formed in the outflows of stars by precipitating on pre- 
existing grains and that MgS in all likelyhood forms as mantle 
on the silicon carbide grains formed closer to the star at higher 
temperature. 

The plan of the paper is as follows: First, we discuss some 
aspects of the chemistry associated with sulphur bearing con- 
densed phases in Sect. [2] The kinetics of dust particle growth 
in a stellar outflow is discussed in Sect. [3] The results of some 
model calculations for MgS condensation in stellar winds are 
presented in Sect. [4] and the optical properties of grains formed 
by MgS grown on a silicon carbide or carbon core are briefly 
discussed in Sect. [5] Our final conclusions are presented in Sect, 
i 

2. Stability of sulphur compounds 

We start with a brief discussion of the equilibrium chemistry of 
sulphur compounds in a carbon rich element mixture. Figure [TJ 
shows the results of a calculation of the chemical equilibrium 
composition of an element mixture with cosmic element abun- 
dances (e.g. Grevesse & SauvaQ 1998) with carbon abundance 
set to 1.2 times the oxygen abundance. The lines in this P-T- 
diagram correspond to the upper stability limits of the indicated 
solid phases. For comparison we also show the wind trajectory 
of a stationary dust-driven wind model for the outflow from a 
carbon star. In unconditional chemical equilibrium also some 
oxygen bearing condensed phases would exist at temperatures 
below 800 K (cf. Lattimer et al. |19781 l. The oxygen for their for- 
mation has to come from a reduction of CO. It is unlikely that 
such reactions are kinetically possible and we assumed here that 
the oxygen is completely blocked in CO which does not react 
with other species in the system, i.e., we consider a conditional 
chemical equilibrium where O is arr ested in CCS 

1 We prefer these over the more recent data of Asplund et al. (2005 ) 
since the reduction of the C and O abundances in the new table are in 
clear contradiction to the results of helioseismology (e.g. Chaplin et al. 

[2007T > 

2 Some condensation temperatures then are somewhat different from 
the results of more conventional condensation calculations 



From the abundant elements with abundances relative to H 
exceeding 10~ 6 , only CaS and MgS are formed as solid sulphur 
bearing phases in chemical equilibrium. Troilite (FeS), that is 
abundant in oxygen rich environments, is not formed in the car- 
bon rich environment because MgS is more stable than FeS, and 
in a carbon rich environment the Mg is not consumed by the 
even more stable magnesium silicates as in the oxygen rich case. 
The first sulphur bearing solid compound that appears in a cool- 
ing sequence at pressure 10~ 10 bar is CaS at about 950 K, fol- 
lowed by MgS at about 750 K. The temperatures where CaS and 
MgS become stable in a stellar outflow depend somewhat on the 
mass-loss rate and are of the just mentioned order of magnitude 
for a mass-loss rate of M — 10~ 5 M0 yr -1 or somewhat less for 
smaller values of M. 

Calcium sulphide may be formed in a carbon rich environ- 
ment prior to MgS condensation, but because of the low calcium 
abundance (which is only 1/17 of the Mg abundance) it would 
be difficult to detect this dust component in infrared spectra from 
circumstellar dust shells in the presence of MgS dust. The for- 
mation of CaS cannot consume all the S since the sulphur abun- 
dance is about eight times the Ca abundance, i.e., MgS is formed 
additionally to CaS at lower temperatures. Since the Mg abun- 
dance is about twice that of the S abundance, at most half of the 
available Mg can be consumed by the formation of MgS. The 
remaining fraction of Mg stays in the gas phase and may or may 
not precipitate later on other condensed phases. 

The abundance of S bearing molecules in the gas phase is 
shown in Fig.[2]for a cooling sequence with total pressure 10~ 10 
bar that is representative for the pressure in the dust formation 
zone of a stellar wind (cf. Fig.[TJ. The figure shows the composi- 
tion of the gas phase in a restricted chemical equilibrium where 
carbon is not completely condensed into solid carbon, a situa- 
tion encountered in an outflow, where carbon formation cannot 
keep pace with cooling and only part of the carbon is converted 
from gaseous species into the condensed phase. Results for the 
sulphur chemistry are shown for the two cases that (i) the Si con- 
denses into silicon carbide (SiC) at temperatures below about 
1250 K and (ii) for the case that formation of SiC is suppressed. 
A comparison of the results for both cases shows: 
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- If no SiC is formed, the sulphur is almost completely bound 
in SiS below 1500 K and no other S bearing species are abun- 
dant in the gas phase. 

- Only if the Si is consumed by the formation of solid SiC, the 
sulphur is liberated and forms increasing amounts of other 
S-bearing molecules - in particular H2S - in pace with con- 
version of Si from SiS into silicon carbide. 

This means that the sulphur chemistry of the gas phase and the 
formation of S-baring solid phases are strongly coupled to the 
process of silicon carbide condensation. 

As long as not much solid silicon carbide has formed the 
sulphur is blocked in the rather stable SiS molecules and the 
formation of MgS from abundant gas phase species requires the 
following reaction 

SiS + Mg -^MgS(s) + Si. (1) 

The change of free enthalpy by this reaction is AGi = 
-30kJ/mol at 800 K (thermochemical data are from Chase 
1998 ). After the consumption of SiS by silicon carbide forma- 
tion and formation of abundant H2S the following reaction be- 
comes possible 

H 2 S + Mg MgS(s) + H 2 . (2) 

The change of free enthalpy in this case is AG2 = -350kJ/mol. 
The corresponding laws of mass action for these two chemical 
reactions are 

aMgswPsi = PsisPM g e~ AG,/RT (3) 
aM g s( S )PH 2 = Ph 2 s Pm s e~ AG2/RT (4) 

where flMgS(s) is the activity of solid MgS. In chemical equi- 
librium the activity of the solid compound equals unity. These 
equations determine the equlibrium pressure in Sect. 13.21 

3. Formation of MgS in a stellar wind 

In stellar outflows temperature and density variations occur on 
shorter timescales than particle growth. In this case condensa- 
tion proceeds far from chemical equilibrium and particle growth 
has to be treated by reaction kinetics. We consider here the con- 
ditions under which MgS condensation may occur in a rapidly 
expanding stellar outflow. 

3.1. The wind model 

3.1.1. Stationary wind model 

The formation of MgS in stellar outflows from AGB stars seems 
to be associated with late phases of AGB evolution of low and in- 
termediate mass stars, where the mass-loss rate is very high and 
the stars have lost already most of their hydrogen rich envelope. 

At the tip of the AGB, radiation pressure on dust be- 
comes the dominating driving source of the wind. The un- 
derlying stars seem all to be variables, either Miras or LPVs. 
The structure of the outflow is complicated in this case, since 
there are always several shocks running outwards, which are 
superposed on the average outflow of stellar material. This is 
demonstrated by several model calculations of dust forming stel- 
lar outflows of pulsating stars (e.g. Bowen I1988I Fleischer et 
al. [1991] [T9921 Fe uchtinger et a l. IT9931 Hof ner & Dor fi [19971 
Hofner et al. [19981 Winters et al. 1997 [20001 Jeong et al. |2003l . 



Calculating dust formation rates for multi-component dust 
mixtures of stars with quite different elemental compositions 
and widely varying stellar parameters presently seems not to 
be possible for reasons of computational time requirements if 
the shock structure of the wind is to be taken into account. For 
the purpose of calculating the composition of the dust mixture 
and the amount of dust formed in the outflows of AGB-stars we 
use a more simple model, which assumes a stationary outflow. If 
one compares the velocity and density profiles of such station- 
ary winds with published models of dust forming pulsators (e.g. 
Winters et al. 2000, Jeong et al. 2003 ), one observes in most 
cases a strong resemblance of the average velocity and density 
profiles with that of the stationary models, except that in the pul- 
sation models outwards propagating shocks are superposed on 
the average flow structure. So we can hope to obtain at least an 
estimate of the average quantities of dust formed in the outflow 
if dust formation is calculated for such an average outflow struc- 
ture. One important consequence of this neglect of the detailed 
structure of the flow is that we cannot determine self-consistently 
the mass-loss rate M. We have to treat this as a free parameter, 
which has to be fixed in some other way. 

3.1.2. Stellar parameters 

The central star is assumed to be on the tip of the AGB. Its lumi- 
nosity then is typically L» = 1.5 X 10 4 Lq. 

For stars forming MgS in their outflow a rather high metal- 
licity is required since at low metallicities the abundances of Mg 
and S would be too low for this to be possible. Only a small 
fraction of the carbon stars in the Small Magellanic Cloud, for 
instance, seem to form MgS dust (Sloan et al. 2006 ), while for 
carbon stars in the Large Magellanic Cloud the formation of 
MgS seems to be more common (Zijlstra et al. 2006 ), and in the 
Milky Way it is quite common for highly evolved carbon stars 
(e.g. Hony |2002l . We note already at this place that according 
to observations MgS formation seems to be associated with SiC 
formation: Stars showing the MgS feature also show the SiC fea- 
ture, but not all stars showing a SiC feature also show the MgS 
feature in their spectra (e.g. Zijlstra 2006 ). In order that low and 
intermediate mass stars of Pop I (or not much lower) metallicity 
become carbon stars, their initial masses must be from the mass 
range between about 1.5 Mq and about 4 Mq. 

After massive mass-loss on the AGB the stellar mass is likely 
to be already substantially reduced and we assume in the model 
calculation a fixed stellar mass of 1 M0. The stellar radiation is 
approximated by a black body spectrum. The effective tempera- 
ture (in the inter-pulse phase) is determined from the relation 

log r eff = .234 log M - .2 log L - .1 16 log ^ + 4.146 (5) 

given by Vassiliadis & Wood d 19931 ) for stars on the AGB. The 
mass M and luminosity L are in solar units. With the assumed 
stellar parameters and solar metallicity one finds T e s = 2 250 K. 
The mass-loss rate is estimated from the relation 

M = 2.1xl0~ 8 — (6) 

Voo 

given by Vassiliadis & Wood ( 11993b for stars in the superwind 
phase. The luminosity is in units of Lq, the outflow velocity Voo 
in kms -1 , and the mass-loss rate in units of M0yr~'. With the 
estimated stellar parameters and a typical outflow velocity in the 
superwind phase of 10 . . . 15 km s one obtains a mass-loss rate 
of about M = 2 . . . 3 x 10~ 5 M© yr _1 for stars at the end of their 
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AGB evolution. A value of M = 3 x 10~ 5 Mq yr~' is used in the Table 1. Basic data used for calculation of MgS grain growth 
model calculations. 

A typical run of density and temperature in a stationary stel- 
lar wind model of a carbon star calculated for a dust driven sta- 
tionary wind model (cf. Gail & Sedlmayr [19861 [19871 for the 
details of how such models are constructed) is shown in Fig.Q] 

The C/O abundance ratio of the outflowing material in any 
case exceeds unity since MgS is not formed in an oxygen rich 
environment. 

3.2. Dust particle growth 

The model calculation of condensation in the outflow from a car- 
bon star considers the following dust components: Carbon dust, 
silicon carbide dust, and iron dust. The formation of these dust 
species is calculated as outlined in Ferrarotti & Gail (2006). Here 
we add the growth of magnesium sulphide to our model. The 
dust formation is calculated by solving the equations for dust 
growth (e.g. Gail 120031 
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for all species of interest. Here a is the radius of a dust grain, Vq 
the volume of one formula uni{3 of the chemical compound in 
the solid, a is the growth coefficient, n gr is the particle density of 
the growth species in the gas phase, m gr its mass, p eq is the partial 
pressure of the growth species in a state of chemical equilibrium 
between the condensed phase and the gas phase, and T g and 7d 
are the gas temperature and the internal lattice temperature of 
the dust, respectively. The quantity 



(D = 



np.m H 
1 + w z 



%kT s 



(8) 



depends in principle on the drift velocity w of the grains relative 
to the gas. Here drift is neglected and we put <D = 1 since for high 
mass-loss rates drift velocities are small compared to the sonic 
velocity. The volume Vq of one chemical formula unit occupied 
in the solid is 



V 



Qd 



(9) 



where is the atomic weight of the chemical formula unit and 
Qi the bulk density of the condensate. 

The growth equation has to be complemented by an equation 
for the consumption of the growth species from the gas phase. 
The details are described in the paper by Ferrarotti & Gail (2006) 
and are not repeated here. 

As rate determining step for the formation of MgS via reac- 
tion © we assume the addition of a H2S molecule from the gas 
phase since this is the least abundant of the molecules involved 
in the reaction. p eq in Eq. (0 is calculated from from Eq. (0]l 
with activity «MgS = 1 ■ The thermodynamic data for calculating 
AG are taken from JANAF-tables (Chase [19981 . 

The growth coefficient a, the probability that a collision is 
followed by a growth step, seems not to be known. For the simi- 
lar solid MgO with the same structure and similar bonding prop- 
erties Hashimoto ( 1990) measured a growth coefficient of about 
0.2 and we will use the same value for MgS. 



3 The formula unit corresponds to the chemical formula of the con- 
densed phase, if it is written with integer stoichiometric coefficients 



3.3. Conditions for grain growth 

Some insight into the origin of the difficulty to explain how MgS 
forms in the outflow can be gained from considering a simplified 
growth model. We neglect in the growth equation (i) the vapor- 
isation and (ii) the consumption of the growth species by grain 
growth. Additionally we use in the expression for the thermal 
velocity a representative value for the temperature in the growth 
zone of the dust. The first neglect is generally not too bad for dust 
growth problems since the vapour pressure p eq usually rapidly 
decreases with decreasing temperature below the stability limit 
of a solid. The second neglect restricts the applicability of our 
considerations to the early growth phase before the gas phase 
becomes exhausted of condensable material. This initial growth 
phase is that which we now consider. 

The particle density n gl is the least abundant gas phase 
species involved in the rate determining reaction for the growth. 
Its abundance is usually some fraction g of the least abundant el- 
ement from which this gas phase species is composed of. In our 
case, where reaction (f2]) is assumed to be the key reaction for 
MgS growth, the least abundant growth species is H2S and sul- 
phur is the least abundant of the elements forming this molecule. 
We write « gr = ge e \NH, where A^h is the number density of H nu- 
clei which can be calculated from the mass density p of the wind 



P 



(1 +4e He )m H 



(10) 



The mass density is given in a stationary outflow with mass-loss 
rate M and outflow velocity v by 



P = 



M 
4nr 2 v 



(11) 



Since condensation of MgS becomes possible only at a tem- 
perature much lower than the condensation temperature of car- 
bon dust, the radiation pressure on carbon dust has driven the 
outflow velocity already nearly to the terminal outflow veloc- 
ity of the stellar wind if MgS condensation commences. For the 
problem of MgS condensation we can consider the wind velocity 
v therefore to be constant. By changing the independent variable 
t to r by means of vdt = dr we can write the equation for the 
initial growth of dust as 



da 

— = Vq a Vgr 

dr 



M 



1 Amu 47zt 2 v 2 



(12) 
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Here v gr is the thermal velocity of the growth species (assumed 3.4. Heterogeneous growth 



to be constant). Integrating from the inner radius /?d of the onset 
of growth of the considered dust species to infinity we obtain 



a + V a Vgr ge e \ 



M 



1 .4m H 4nRa v- 



2 ' 



(13) 



where ciq is the radius of the seed nuclei for dust growth and a M 
the radius at infinity. This equation requires that the grain radius 
floo remains smaller than the maximum radius a max attained if 
all condensable material is condensed. This maximum possible 
radius is given by 



Vofkey 
Vkey 



(14) 



Here ed is the number of dust grains per hydrogen nucleus, ek ey 
the key element for the formation of the dust species considered, 
and Vkey the number of atoms of the key element in the chemical 
formula of the solid. 
We can now write 



O.Q + a n 



M 



with 



flmai 1 -4m H 4nR A V 2 
VoffVgrgeel 



(15) 



(16) 



The quantity M cr defines a critical mass-loss rate for particle 
growth. If the mass-loss rate M of the stellar outflow is lower 
than M cr , the final particle radius of dust grains of the considered 
kind in the outflow is smaller than the maximum radius attained 
in case of complete condensation. Condensation by growth start- 
ing with some kind of small seed particles is always incomplete 
in this case and most of the condensable material remains in the 
gas phase. If, however, the mass-loss rate M of the stellar out- 
flow is of the order of M cl - or exceeds M CI , the final particle radius 
in the outflow approaches a max ; complete condensation is possi- 
ble in this case by growth starting with some kind of small seed 
particles (at least this is not kinetically forbidden). 
Applying this to MgS condensation we obtain 



cr.MgS 



5.4 x 1(T 4 Moyi-- 1 



(17) 



Here we assumed a typical outflow velocity of v = 10 kms -1 
and a radius of /?d = 2 x 10 14 cm corresponding to about 5 stellar 
radii. As key element for MgS condensation we assume sulphur 
since it has lower abundance than Mg. The growth species is 
H2S. As typical temperature for calculating v gr we choose 600 K. 
All other quantities are given in Table Q] The critical mass-loss 
rate is about a factor of ten higher than the highest mass-loss 
rates observed for highly evolved AGB stars (e.g. van Loon et al. 
1999). As a result the radius can approach at most 1/10 of that 
for complete condensation and only about 1/1000 of the con- 
densable material can be condensed into MgS grains if the dust 
grains start their growth for instance from seed particles with 
sizes of the order of 1 nm that are typical sizes for seed particles 
resulting from nucleation from the gas phase. 

For such a low fraction of condensation of S and Mg into 
MgS this dust component would not be detectable in the infrared 
emission from circumstellar dust shells of AGB stars. Since MgS 
is frequently observed in highly evolved AGB stars, it cannot be 
formed by condensation from the gas phase via nucleation and 
subsequent growth. We need therefore to look for some other 
kind of formation mechanism. 



The most promising alternative is heterogeneous growth on 
some other kind of dust particles that have already not-too-small 
radii at the onset of MgS growth in order that the frequency of 
collisions of such a particle with growth species from the gas 
phase is high already from the beginning on. We derive the con- 
dition for this type of growth to occur by considering the initial 
growth phase, where the growth of a coating of MgS on some 
type of carrier grain has not yet significantly increased the total 
grain radius. 

We use essentially the same type of simplifications as before 
and consider the following equation for the increase of the vol- 
ume V of the coating on the surface of carrier grains of radius a c 



dV 2 
— = V av gr n gT 4na c 



(18) 



Introducing again r as independent variable we obtain upon inte- 
grating this equation from the radius R& of the onset of formation 
of the coating to infinity 



M 



Voc = V av gI 4na- c ge A — 

1 Amu 4nRd v L 



(19) 



This is the final volume of the coating at large distances. The 
maximum possible volume V max of the coating on a carrier dust 
grain is, again, given by Eq. ( TBI . We can now write 



Voo = Vrr 



M 

M het 



where 



M het = 



ekey lAmuRdv 



Vkey^d CVgr Ac £ e el 



(20) 



(21) 



is a characteristic mass-loss rate of the stellar outflow for forma- 
tion of coatings on the surface of already existing grains. If the 
mass-loss rate M of the stellar outflow is less than the charac- 
teristic mass-loss rate Mhet only part of the condensable material 
will condense in form of a coating on the surface on some car- 
rier grains. If the mass-loss rate M of the stellar outflow exceeds 
the characteristic mass-loss rate Mhet the volume of the coating 
approaches the maximum possible volume and complete con- 
densation is possible, at least in principle. 

Inserting numerical values for MgS condensation and as- 
suming that the carrier grains for the formation of the coating 
have a radius of a c = 0.1 fim one obtains 



M het = 2.8 x lO" 4 I Moyr- 1 



0.1 /an \ 
a c J 



(22) 



This characteristic mass-loss rate, again, is about a factor of ten 
higher than the highest mass-loss rates observed for AGB stars, 
but since it now refers to the volume we have the result, that 
at the end of the AGB evolution up to about 10% of the con- 
densable material for MgS formation can really be condensed as 
coating on suited carrier grains, if such exist. If the core grains 
have a bigger radius than the assumed 0. 1 /mi, or if the mass- 
loss rate should be exceptionally high the fraction of condensed 
material may be even bigger. Such amounts of condensed MgS 
should be clearly visible in the infrared spectrum. 

Hence, we come to the conclusion that the observed MgS in 
the outflows from carbon stars on the AGB does not form as a 
separate dust species but condenses as coating on the surface on 
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Fig. 3. Evolution of MgS and CaS grain radius as function of 
distance from the centre of the star, for MgS formation (i) on 
seed nuclei of 1 nm size, and (ii) as mantle on a CaS core. Also 
shown are the results of a model with delayed onset of MgS 
condensation on CaS core; details are given in the text. 
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Fig. 4. Evolution of radius of SiC and MgS grains as predicted 
(i) by the model for MgS condensation on seed nuclei, and (ii) 
as mantle on a SiC core (thin and thick lines, respectively). 



some different kind of carrier grains. The formation of signifi- 
cant amounts of MgS dust requires mass-loss rates of at least a 
few times 10~ 5 Mq yr _1 , i.e., this is possible only during the su- 
perwind phase short before termination of AGB evolution if the 
mass-loss rates have grown to very high values. This is just what 
one seems to observe. 

As a side-remark we note the following: Since the available 
amount of S does not suffice for the formation of thick coatings 
of MgS such that the final volume is much bigger than the vol- 
ume of the carrier grains, the volume of condensed material on 
the surface of the carrier grains increases quadratically with the 
radius of the carrier grains (cf. Eq. |19t . Therefore it may well 
be that most of the MgS is concentrated on a sub-population of 
lager-than-the-average grains, if one has a broad distribution of 
grain sizes for the possible carrier grains. 



3.5. Possible carrier grains 

The dust species that may serve as carrier grains for MgS forma- 
tion need to form dust grains with a size not much smaller than 



about a c = 0.1 /mi because otherwise the condensed fraction of 
MgS becomes too small if the radius of the cores is smaller than 
this. This restricts the possibilities, which dust species may serve 
as substrates of MgS growth, to the most abundant dust species. 

The outflowing gas enters the temperature regime below 
about 650 K favourable for MgS condensation with two major 
dust species that have grown farther inside in the shell to the re- 
quired size: SiC grains and carbon grains. From investigations of 
presolar dust grains it is known that from the many poly-types of 
SiC only two are formed in circumstellar environments (Daulton 
et al. 120031 1: The cubic 3C poly-type (usually denoted as /?-SiC) 
is found in 79.4% of all observed grains, and the hexagonal 2H 
poly-type (usually denoted as a-SiC) is found in 2.7 % of all 
observed grains. In 17.1 % of all grains the material is an inter- 
growth of these two poly-types. 

The lattice structure of the SiC grains in stellar outflows, 
thus, in most cases is cubic. Magnesium sulphide also has cu- 
bic lattice structure. The electronegativity difference (cf. Pauling 
|T960j > of MgS is 1.3 and that of SiC is 0.7. Both compounds 
therefore show a significant ionic contribution to their bonding. 
There are obvious similarities of structure and bonding proper- 
ties of SiC and MgS. On the other hand, no such similarity exists 
between MgS and solid carbon. This suggests that SiC is much 
better suited as substrate for MgS precipitation than carbon. We 
propose therefore that MgS grows on SiC grains. 

A possible core grain is any compound that has a lat- 
tice structure similar to the MgS crystalline structure and a 
higher condensation temperature. Among the compounds pos- 
sibly formed in the stellar outflows also CaS satisfies both con- 
ditions. It has a crystal structure, that is the closest to that of 
MgS and becomes stable at a more than 200 K higher tempera- 
ture than MgS (Fig. [T). The abundance of Ca in a stellar wind, 
however, is 16 times lower than the Mg abundance, which is 
not very favourable for its role as core grain for MgS growth. 
Nevertheless, we will also check this possibility. 



4. Models with core-mantle grains 

We performed model calculations of dust formation in stellar 
outflows according to the wind model described in Sect. [3] to 
check the assumption of the core-mantle growth scenario of MgS 
formation. Grain growth is calculated as described in Sect. 13.21 

First we check the model of MgS condensation on tiny seed 
nuclei. The existence of such seeds of 1 nm is taken for granted 
and the growth of MgS is calculated. This results in grain sizes 
of only about 0.02 /jm and a very low fraction of Mg bound in 
dust, the condensation degree, of about 0.01. This is in line with 
our analytical findings in Sect. 13.31 The evolution of grain radius 
and condensation degree for this kind of model are shown in 
Fig. [3] and Fig. [5] respectively. The small amount of condensed 
MgS found in this case is obviously not enough to explain the 
observed 30 /mi emission in the spectra of carbon stars, under- 
lining, again, the need for a core-mantle growth mechanism of 
MgS formation. 

Next we check an alternative growth scenario for MgS for- 
mation, in which first CaS condenses on seed nuclei by the ana- 
logue of reaction ©, and then MgS grows as a mantle on this 
core. The model calculations show that the size of CaS grains is 
much smaller than the required radius of 0.1 /mi. We also tested 
the case, that MgS condensation is delayed by imposing the ad- 
ditional condition for the growth rate /grjvigs ^ 3/ gr .caS to gain 
additional time for CaS core growth in order to get bigger grains. 
This artificial delay results, indeed, in bigger CaS cores, but also 
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Fig. 5. Evolution of condensation degree of Mg with distance 
from the centre of the star for condensation of MgS (i) as mantle 
on top of a SiC core and (ii) on seed particles (solid and dashed 
line, respectively). 

this does not help to increase the final size of the MgS man- 
tle. The results of these calculations of the evolution of the MgS 
grain radius with distance from the centre of a star for normal 
and delayed condensation on CaS are also shown in Fig. [3] The 
figure shows that the condensation of MgS on a CaS core re- 
sults in almost the same grain radius as the growth of MgS on 
seeds grains. Both these formation mechanisms can not provide 
the amount of MgS dust that is observed to condense in stellar 
outflows. 

SiC grain growth, however, commences at higher tempera- 
ture and proceeds faster than CaS growth, so that at the onset of 
MgS formation the SiC core radius reaches oc 0.07 /mi as shown 
in Fig. [4] The results of MgS condensation on tiny seed nuclei 
are shown in Fig.|4]for comparison. The SiC core is big enough 
to allow the formation of a considerable thickness of the MgS 
coating, of order 0.02 /mi, which results in a final grain size of 
0.9 /mi. This value is very close to the required grain sizes to ex- 
plain the shape of the extinction feature in the spectra of carbon 
stars. We assume in the calculation that SiC grains do not grow 
anymore as soon as the formation of the MgS mantle begins. 
Figure |4] show that growth of SiC grain with and without termi- 
nation due to MgS mantle formation does not change the grain 
size very much. Thus, a model with MgS mantle formation on 
SiC does not change noticeably the fraction of Si condensed in 
SiC dust, but results in a much higher value of the degree of con- 
densation of Mg in MgS as compared to the case of growth on 
seeds. This is illustrated in Fig. [5] 

A similar test calculations for MgS coating formation on 
carbon grains also gave sufficient thicknesses of MgS layers. 
Basically, condensation of MgS on carbon cores cannot be ex- 
cluded, but from a physical point of view, MgS mantle formation 
is more likely on a SiC core, since the bond lengths and bonding 
properties of SiC are similar to that of MgS. Laboratory exper- 
iments would be necessary to arrive at a final conclusion about 
the carriers of MgS mantles. 

5. Optical properties of grains with SiC core and 
MgS mantle 

In order to study the difference of optical properties of grains 
with coatings of MgS either on a SiC core or on a carbon core 




"^abs 



Of 



Fig. 6. Wavelength variation of the absorption efficiency C\ 
spherical grains with SiC core of radius r c coated with MgS man- 
tles of thickness r m for a total grain radius r c + r m of 0.1 /urn and 
for MgS-mantle thickness to total radius ratios from to 1 in 
steps of 0. 1 (from bottom to top). 



we calculate the absorption efficiency CJ) bs (i.e. the ratio of ab- 
sorption to geometrical cross-section) in the small particle limit 
(cf. Bohren & Huffman [19831 



with 



; 4xlm- 



4n(r c + r m ) 3 



(23) 



a = 4n(r c + r m ) 



0c + r m y(e 2 - lXg + 2e 2 ) + r c 3 (2e 2 + - e 2 ) 
(r c + r m )He 2 + 2)(e, + 2e 2 ) + r c 3 (2e 2 - 2)(e, - e 2 ) 



(24) 



Here r c is the core radius, r m the thickness of the mantle and e\ 
and e 2 are the complex dielectric functions of the core and man- 
tle material, respectively, x is the size parameter 27r(r c + r m )/A. 
The small particle approximation x «: 1 is valid in our case 
since we are interested in wavelengths from the region from 10 
to 40 /mi that are much bigger than the size of circumsteilar dust 
particles. 

Optical constants for MgS are taken from Begemann et al. 
(1 1994| l for Mg x Fei_ x S with x = 0.9. For carbon we use data for 
the BE carbon dust of the evaluation of the data of Colangeli 
et al. ( [19931 by Zubko et al (fl996l >. For SiC we used the data 
for SiC from Laor & Draine ( 1993). Extinction efficiencies have 
been calculated for SiC grains coated with MgS mantles and, for 
comparison, also for carbon grains with MgS mantles. 

Figure [6] shows results for grains of total radius r c + r m 
of 0.1 fim and mantle thicknesses between and 0.1 /mi, i.e., 
for grains ranging in composition from pure SiC to pure MgS 
grains. The results show the rather broad MgS absorption band 
centred around about 26 /mi, which in itself shows some struc- 
ture. Additionally there appears a secondary peak centred around 
~ 33 ... 38 /mi for not too big mantle thicknesses. Its position 
depends on the size ratio of core and mantle. This feature is not 
present in the extinction efficiencies of MgS or SiC and results 
from the particular run of the complex dielectric functions of the 
core and the coating in this wavelength regiorQ. 



4 It has been checked if this results from a Frohlich mode in coated 
particles at a frequency where the denominator in Eq. i24\ vanishes, cf. 
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Fig. 7. Wavelength variation of the absorption efficiency C ^ of 
spherical grains with SiC (top) or carbon (bottom) core and MgS 
mantle with total radius of 0.1 /mi. The MgS mantles range in 
thickness from zero to unity in steps of 0.05 for the ratio of man- 
tle thickness r m to total radius r c + r m , where r m is the core radius. 
The extinction properties vary from that of pure SiC or carbon 
grains to that of pure MgS grains. An extra extinction feature 
centred on 33 //m for MgS grains with SiC core shows up at 
moderate mantle thickness. This is missing in case of carbon 
cores. 



Figure|7]shows the variation of the extinction coefficient with 
varying core to mantle size ratios - both for MgS grains with 
SiC and with carbon core - in a different representation and an 
extended wavelength region that also covers the 1 1 /mi feature 
of SiC. The total grain radius r c + r m , again, is 0.1 /mi and the 
fraction r m /(r c + r m ) is varied between and 1. In the upper part 
of the picture one recognises how the SiC feature disappears and 
the MgS feature appears as the composition of the grain varies 
from pure SiC to pure MgS. The extra feature peaking at about 
33 . . . 38 /mi, depending on the mantle thickness, is clearly seen 
in case of coatings of moderate, but not too small fraction of the 
total size. For thick coatings the extra feature disappears together 
with the 1 1 /mi feature from the SiC core. The feature is not very 
distinct and will probably hardly be detectable in many cases. 

The particular feature peaking at about 33 . . . 38 /mi is miss- 
ing in the extinction of carbon grains coated with MgS, as can 
be seen from the lower part of Fig. [7] because of the completely 
different properties of the dielectric function. Therefore this ex- 
tra feature can be considered as an indicator of a silicon carbide 



the discussion of this effect in Bohren & Huffman ( 1983 1, but this seems 
not to be the case 



core of MgS grains, if the presence of MgS is indicated by its 
absorption band around 26 //m. Observationally it was found by 
Volk et al. ( 120021 1 that the so called '30/mi feature' of carbon 
stars shows in some cases some structure and seems to consist 
of two overlapping features at 26 /mi and 33/mi. Our results for 
the extinction of SiC-core-MgS-mantle grains indicate, that one 
just sees that the MgS in outflows from carbon stars grows as 
mantle on silicon carbide cores. If it is not seen this probably 
does not mean that in this case MgS precipitates on carbon, but 
merely that the feature is not sufficiently well defined to be un- 
ambiguously be detectable. 

6. Concluding remarks 

In this paper the mechanism of MgS formation in stellar outflows 
of AGB stars is studied. From some elementary considerations 
on the kinetics of MgS condensation we estimate a critical mass- 
loss rate of about 5 x 10~ 4 M0 yr~' for efficient MgS formation 
if this would condense as a separate dust species via nucleation 
and subsequent growth. This value of the mass-loss rate is much 
higher than any observed mass-loss rate of AGB stars. Within 
the observed range of mass-loss rates of AGB stars, however, 
the amount of MgS that can condense on seed particles in stellar 
winds is much too low to explain the observed strong emission 
band from MgS in the infrared spectra of many carbon stars. 

A model of heterogeneous growth of MgS on SiC precur- 
sor grains is proposed here as the most promising mechanism 
of MgS formation. Analysis of the equilibrium chemistry of sul- 
phur compounds in a carbon rich mixture shows that the forma- 
tion of MgS is strongly coupled to the process of silicon carbide 
condensation. If no SiC is formed in a stellar outflow, the sulphur 
is almost completely bound in SiS molecules. Only after signif- 
icant consumption of the abundant gas phase SiS by SiC forma- 
tion, H2S is formed in abundance and reactions of this with Mg 
makes possible MgS condensation. 

Our model calculations of MgS formation in stellar outflows 
show that only MgS growth as mantle atop SiC cores, that con- 
densed before the on-set of MgS formation, results in sufficiently 
high degrees of condensation of Mg into MgS of the order of 
~ 0.2 which are required that MgS may be clearly visible in the 
infrared spectrum. Condensation of MgS on tiny seed particles 
would result in a fraction of the Mg condensed in MgS dust of 
the order of only 0.01. This is much too small to be detectable. 

Additionally, we performed calculations of the extinction 
properties of grains with SiC core and MgS coating of various 
thicknesses. The presence of a SiC core inside of MgS grains 
results within some range of core/mantle volume ratios in a sec- 
ondary peak near 33 . . . 35 /mi in the broad emission band as- 
sociated with MgS, that seems to be observed in some spectra 
of AGB stars. This feature is absent if MgS forms a mantle on 
carbon grains. We propose this feature as an indicator for the 
presence of a silicon carbide core of MgS grains, if the presence 
of MgS is indicated by its absorption band around 26 /mi. 
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